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The Orion spacecraft will be required to perform a three-burn trans-Earth injec- 
tion (TEI) maneuver sequence to return to Earth from low lunar orbit (LLO). The 
origin of this approach lies in the Constellation Program requirements for access 
to any lunar landing site location combined with anytime lunar departure (which 
could require up to a 90° plane change). This paper documents the development 
of optimized trajectory databases used to rapidly model the performance require- 
ments of the TEI three-bum sequence for an extremely large number of mission 
cases. It also discusses performance results for lunar departures covering a com- 
plete 18.6 year lunar nodal cycle as well as general characteristics of the optimized 
three-bum TEI sequence. 


INTRODUCTION 

Use of a three-burn maneuver sequence to transfer from a circular orbit to an escape hyperbola 
has been studied since 1959 when it was first proposed as a method to save propellant over one 
and two-burn escape sequences. 1 During the Apollo lunar program, this approach was investigated 
as a possible maneuver sequence to use for trans-Earth injection (TEI), 2-10 and the reverse prob- 
lem (transfer from a hyperbolic to a circular orbit) was also studied for lunar orbit insertion (LOI) 
applications. 11,12 The three-burn escape maneuver sequence was also studied for applications to 
interplanetary missions. 13 These early studies provided valuable insight into the problem, but made 
many simplifying assumptions in order to obtain performance data in a reasonable amount of time. 
The most common assumption was to ignore the influence of the Earth’s gravity during the maneu- 
ver sequence. Ultimately, the three-burn sequence was not employed during the Apollo program. 

The three-burn TEI problem received renewed interest in 2004 with the inception of the NASA 
Constellation program to return humans to the Moon. 14-17 The current Constellation lunar architec- 
ture requires the Orion vehicle to provide access to the entire lunar surface with an anytime abort 
return capability. This requirement means that the vehicle must be capable of executing up to a 90° 
plane change to transfer from low lunar orbit (LLO) to the hyperbolic Earth return trajectory. In 
order to make this possible without a prohibitively large propellant cost, a three-burn TEI maneuver 
sequence is necessary. Initial Constellation work employed many of the same simplifications used 
in the Apollo-era studies, but as the Orion vehicle design matured, more refined estimates of flight 
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performance became essential since they were being used to size the propellant tanks. Because the 
Constellation architecture is intended to be operational for a number of years, the entire range of 
possible Earth-Moon geometries must also be considered. A statistical approach to performance as- 
sessment was also desired, so that the vehicles are not over-designed to account for an insignificantly 
small percentage of worst-case mission scenarios. 

At the NASA Johnson Space Center (JSC) Flight Mechanics and Trajectory Design Branch, the 
Copernicus trajectory design and optimization system 18 has been used extensively to study the TEI 
sequence for Constellation and Orion work*. Copernicus uses a multiple shooting parameter opti- 
mization approach with explicitly-integrated trajectory segments. It can be used to optimize end-to- 
end trajectories from low Earth orbit (LEO) to LLO, and from LLO to entry interface at the Earth, 
including complex environmental effects such as multi-body perturbations, high order Earth and 
Moon gravity models, and finite burn maneuvers. 

The convergence time for fully optimized TEI trajectories represents a limiting factor on the num- 
ber of cases which can practically be run. Given the range of mission variables (e.g., flight times, 
departure orbits, and entry constraints), it is not yet practical to generate fully optimized high- 
resolution end-to-end trajectories for billions of cases across an entire lunar nodal cycle. Therefore, 
a database approach was developed. 16 Strategic runs exploiting the high resolution trajectory op- 
timization capability of Copernicus were used to construct databases of the performance costs for 
specific mission phases (such as the three-burn TEI sequence). These databases could then be used 
to quickly “look up” the performance of the various mission phases for specific cases, allowing the 
assessment of a huge number of mission permutations very quickly. 

BACKGROUND AND FUNDAMENTALS 
The Three-Burn Sequence 

The two fundamental principles behind the three-burn TEI escape sequence are: (1) the direction 
of motion is most efficiently changed where the velocity is lowest, and (2) the energy level is most 
efficiently changed where the velocity is highest. 2 A schematic of the maneuver sequence is shown 
in Figure 1(a). The first maneuver (TEI- 1) takes place in the initial circular orbit (currently baselined 
for Orion at 100 km altitude) and inserts the spacecraft into a large elliptical orbit about the Moon. 
The second maneuver (TEI-2) occurs near apoapsis of this intermediate orbit and performs most of 
the plane change. The third maneuver (TEI-3) takes place near periapsis of the second intermediate 
orbit and inserts the spacecraft onto the departure trajectory that returns to Earth. The TEI-1 and 
TEI-3 maneuvers are mostly coplanar, although they can each perform a small amount of plane 
change. The departure trajectory is hyperbolic with respect to the Moon, and elliptical with respect 
to the Earth. For Orion, the pre-TEI-3 periapsis altitude is limited to at least 100 km as a “fail safe” 
factor, to ensure that a failed TEI-3 will not result in a surface impact (since the optimal TEI-3 will 
generally occur just before periapsis of the second intermediate ellipse). 

For Orion mission trade studies, the total flight time from TEI-1 to TEI-3 (A txEi) has been 
limited to a maximum of 48 hours. With the current vehicle design, the most stressing missions will 
require a 48 hr TEI duration from TEI-1 to TEI-3, allowing the plane change maneuver to occur as 
far from the Moon as possible. For easier missions, where propellant margin is available, shorter 
flight times can be used. Operationally, other time constraints may be necessary to account for 
other factors, such as crew sleep cycles and navigation updates. When a large plane change is not 

* Recent research involves generation of analytical initial guess trajectories to be used for optimizer initialization 19,20 
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Figure 1 Problem Setup. The full three-burn TEI to El problem is subdivided into 
two separate problems: estimation of the vector using a single-burn departure, 
and a three burn TEI sequence to target the Voo vector. 


Table 1. Mission Cases 


Mission Case 

TEI-1 to TEI-3 (hrs) 

TEI-3 to El (hrs) 

Earth Entry Constraints 

Fast expedited A 

24 

89 

h, 7 

Fast expedited B 

36 

77 

h, 7 

Relaxed expedited 

48 

89 

h, 7 

Simplified Nominal 

48 

89 

h, 7, Az 

Entry Target Nominal 

48 

89-114 

h, 7, Az, (/>, A 



-240 -220 -200 -180 -160 -140 -120 -100 -80 

Longitude (deg) 


Figure 2 Entry Target Line Constraints. Sixth-order polynomials for longitude and 
azimuth as functions of latitude are used to impose entry constraints for landing off 
the coast of San Clemente Island. The target line is an approximation to the entry 
interface points. 
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required, a one-burn or two-burn sequence may be employed. These issues are not considered in 
this study. 

Mission Cases 

Five classes of missions, each defined by different flight times and Earth entry conditions, will 
be considered in this study. These mission cases (summarized in Table 1) are currently being used 
to determine the Orion vehicle tank size requirements. The Earth entry interface (El) state at tEi is 
parameterized as: altitude (h), latitude (</>), longitude (A), velocity (v), azimuth (Az), and flight path 
angle (7). For all cases, the El state component h is specified to be 121.92 km and 7 is specified to 
be -5.86°. 

“Expedited return” missions are defined by El constraints on altitude and flight path angle only, 
with no constraints on the location or direction of the entry point (latitude, longitude, or azimuth). 
They are designed for abort situations (e.g., vehicle or equipment failure, crew injury, etc.) where 
the vehicle must return to Earth quickly and/or compute maneuvers without contact with Earth, 
necessitating the use of simplified on-board non-optimal targeting methods. The fast expedited 
returns allow for a total transfer time (from TEI-1 to El) of 1 13 hrs. There are two cases considered 
here (A and B), which are distinguished by the relative amount of time spent from TEI-1 to TEI-3, 
and from TEI-3 to EL Relaxed expedited returns allow for a total transfer time of 137 hrs. 

“Nominal” missions are more constrained at El in order to target specific landing sites or re- 
gions. 21 For trade studies, the simplified nominal case targets a zero azimuth entry with a total 
transfer time of 137 hrs. This is an approximation of a worst-case targeted entry condition. The 
entry target nominal case enforces a specific relationship between latitude, longitude, and azimuth 
at El in order to target a specific landing site. In this study, the entry target constraints are imposed 
by using two sixth-order polynomials expressing longitude and azimuth as functions of latitude, as 
shown in Figure 2. The polynomials are constructed from data corresponding to valid entry states 
for a U.S. coastal water landing near San Clemente Island. 21 The smooth polynomial constraints 
(rather than the discrete points of the entry map) are more suitable for the optimization process, and 
closely reflect the actual performance requirements. This case requires a variable flight time (from 
137 to 162 hrs) to achieve the specified longitude. 

Lunar Cycle 

The 18.6 year lunar nodal cycle is the period in which the Moon’s line of nodes (the intersection 
of the lunar orbit and ecliptic planes) completes one revolution. The Earth-Moon geometry is a 
critical driver for TEI performance and thus cannot be ignored when considering the performance 
cost of the three-burn sequence. The lunar nodal cycle used in this study is the period from January 
1, 2018 to August 7, 2036. The full cycle allows a complete range of Earth-Moon- Sun geometries 
to be considered across the expected operational lifetime of the Orion vehicle. 

GLOBAL PERFORMANCE DATABASES 
Overview 

The fundamental assumption of this study is that, for performance purposes, the three-burn TEI 
transfer problem can be divided into two sub-problems: (1) computation of a departure Voo vector 
at the Moon corresponding to the desired return trajectory, and (2) computation of a three-burn TEI 
sequence that achieves the desired vector. Parameterized databases for each sub-problem can 
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be constructed and queried for any mission case to produce a TEI Av value which agrees closely 
with the fully-optimized TEI to El sequence. For this study, the databases are generated with the 
Copernicus trajectory optimization system, using the SNOPT optimizer . 22 


Table 2. Fixed Azimuth Voo Database Parameters. 


Parameter 

Range 

Departure epoch 
Az 
At# 

From 2018-01-01 to 2036-12-31 in 12 hr steps 
From -10° to 190° in 5° steps 
From 65 hrs to 1 13 hrs in 12 hr steps 


Table 3. Entry Target Nominal Database Parameters. 


Parameter 

Range 

Departure epoch 

At# 

From 2018-01-01 to 2036-12-31 in 6 hr steps 
2 cases (65-89 hrs, and 89-114 hrs) 


Voo Vector Database 

For a specified departure epoch (£ 3 ), return flight time (At#), and set of El conditions (h, 7, 
and Az), the escape vector can be estimated by optimizing a single-burn TEI departure from 
a circular polar lunar orbit (as shown in Figure 1(b)). The optimization variables for this problem 
are: longitude of the LLO ascending node (£2 llo)> argument of latitude (u) of the maneuver, TEI 
Av magnitude and direction, and the El state parameters 0, A, and v. In the problem setup, the 
entry state is integrated backwards in time (allowing the known El parameters to be specified, rather 
than requiring them to be imposed as constraints), and the post-TEI maneuver state is integrated 
forwards in time. A state (position and velocity) equality constraint is imposed at the point one day 
after TEI to produce a continuous trajectory. The objective function to be minimized is the TEI 
Av magnitude. The entry target case is slightly modified, since At# and Az are also optimization 
variables, and polynomial constraints on </>, A, and Az are imposed on the El state. 

The solution of this optimization problem produces two local minima, corresponding to two 
different Qllo values separated by about 180°. These are referred to as the North and South 
solutions, since the spacecraft departs over either the lunar north or south pole. The vector average 
of these two vectors at a given epoch provides an approximate vector which can be targeted 
from all orbits at that epoch. For the fully optimized TEI to El problem, the vector will vary 
slightly for different departure orbits, shown as an oval- shaped set of points in Figure 3(a), with the 
average vector being a point near the center of this oval. When joined together with the three-burn 
database, the epoch for the single-burn TEI maneuver becomes the epoch of TEI-3, and the 
vector becomes the escape condition after TEI-3. 

Two databases were generated: (1) a 3D database used for the nominal and expedited return 
cases with the following input parameters: TEI-3 departure epoch (£ 3 ), transfer time from TEI-3 to 
El (At#), and Earth entry (El) azimuth, and (2) a ID database for the entry target nominal case, 
where the only input parameter is epoch. The entry target database is generated for specific 24 hr 
return time ranges (in this case, from 65 to 89 hrs, and from 89 to 114 hrs). Each database was 
generated by stepping through the epochs from 2018 to 2036, using the previous epoch as an initial 
guess, solving the optimization problem, and saving the average vector for each case. The 
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run parameters for these database are given in Table 2 and Table 3. Example data from the fixed- 
azimuth database is shown in Figures 4(a) and 4(b). Each case can be solved in about 1 second 
using Copernicus. The fixed azimuth database produced a total of 2,276,484 cases (about 26 days 
of CPU time), and the entry target nominal database produced a total of 1 1 1,044 cases (about 31 hrs 
of CPU time). 


Three-Impulse Departure Vector 
El Azimuth = 0°, TEI 3 Epoch = 03 April 2024 



(a) 


TEI AV Across Voo Relative Declination 


El Azimuth = 0°, TEI 3 Epoch = 03 April 2024 



(b) 


Figure 3 vectors (colors indicate magnitude in m/s) and TEI Av for opti- 

mized three-burn departures to El from a range of orbits for a simplified nominal 
mission case at a single epoch. 


TEI Three-Burn Database with No Third-Body Perturbations 

The starting point for investigating the three-burn TEI problem to target a Voo vector was to 
make the simplifying assumption that only the Moon’s gravity affects the trajectory during the 
maneuver sequence. Under this assumption, the TEI Av can be parameterized using only three 
reference parameters: (1) the magnitude of the vector (VooX (2) the relative declination (a), 
the angle between the Voo vector and the LLO orbit plane, and (3) the time-of-flight from TEI-1 
to TEI-3 (A tTEi)- A database of this type was constructed with Copernicus using the parameter 
ranges given in Table 4 and was used in early Constellation and Orion studies. 16 Given these three 
parameters, the optimization variables for the three-burn TEI problem are: the argument of latitude 
( u ) of TEI-1 in the LLO, the magnitudes and directions of each maneuver (TEI-1, TEI-2, and TEI- 
3), and the intermediate flight time from TEI-1 to TEI-2. A constraint is imposed that the post-TEI-3 
escape hyperbola has the target Voo vector. The objective function to be minimized is the total TEI 
Av magnitude (Aiq + Av 2 + AU3). Given each input parameter permutation, a parking orbit and 
Vqo vector are constructed, and a three-burn TEI sequence is optimized in Copernicus (as shown in 
Figure 1(c)) with only the Moon’s gravity present in the force model. Representative results from 


Table 4. 3D TEI Database Parameters 


Parameter 

Range 

Uoo 

a 

AtTEI 

From 700 to 1500 in 50 m/s steps 
From 0° to 90° in 5° steps 
From 1 day to 2 days in 0.25 day steps 
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Mag (m/s) 


Magnitude of the North Departure Vector 


- Az=0° — * — Az=90° — * — Az=180° 



\ZS/S/S/W v n/x/vWWWSA\/'x 


12/19/2027 02/07/2028 03/28/2028 05/17/2028 07/06/2028 08/25/2028 10/14/2028 12/03/2028 01/22/2029 

Epoch 

(a) Data from the 89 hour return time fixed-azimuth Vx> database shows the variation of the Vx> magnitude over 
a one-year period for different entry azimuths. 

Post-TEI Voo vectors for 89 hr return time 



-400 

400 


y O/s) 


(b) Voo evolution for a one month period for the 89 hour return case (shown in a Moon- fixed reference frame). 


Figure 4 Data from the fixed-azimuth vector database for three entry azimuth 
values (0°, 90°, and 180°). Here, the vectors are expressed in a Moon-fixed frame. 
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Figure 5. Data from the 3D TEI Database shows the TEI Av (m/s) as functions of a, and A tj'Ei* 


this database are shown in Figure 5. Generally, the TEI Av increases with increasing a, increasing 
Voo, or decreasing A trEi- 

For a given mission case and departure epoch, the parameter that has the strongest influence on 
the TEI Av is a , the relative declination angle between the initial orbit and the escape hyperbola. 
Figure 3(b) shows the dependence of TEI Av on relative declination for the full TEI-1 to El problem 
example given in the previous section. It is clear that, while there is a trend based on relative decli- 
nation, other significant factors are also at play. The general effect of Earth gravity perturbations is 
to spread out the data points. The level of uncertainty from neglecting the third-body perturbation 
effect turned out to be up to ±45 m/s for 48 hour TEI flight times, which was unacceptable for 
the purpose of sizing the Orion propellant tanks. However, the original three-parameter database 
captures the general trends of the data very well. It has a very high case density, and is amenable 
to interpolation via splines, so it is very responsive to the curvature and inflections of the data. For 
these reasons, it is still employed as part of the solution, as explained below. 

Design of the TEI Three-Burn Database including Earth Perturbation Effects 

Figure 6 shows the acceleration field due to Earth’s gravity perturbation of an orbit about the 
Moon, expressed as a percentage of the lunar central-force acceleration. The scale is 10% per 
5000 km (1 grid square) on the plot. The Earth is to the left of the plot on the x-axis, 384,400 km 
distant from the Moon. This field is symmetrical about the Earth-Moon axis. Unperturbed one-day 
(green) and two-day (red) orbits are shown for scale, with marks every three hours of travel. The 
perturbation acceleration, expressed as a fraction of the main-body acceleration, grows non-linearly 
with radius. For the two-day orbit, the most significant perturbation effects occur within ±18 hours 
of apoapsis. 
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Figure 6 This plot shows accelerations due to Earth gravity perturbations of an orbit 
about the Moon, expressed as a percent of the central-force accelerations. The scale 
is 10% per 5000 km (1 grid square) on the plot. One day (green) and two day (red) 
orbits are shown for scale. 


In designing the TEI three-burn sequence database that included Earth perturbation effects, the 
challenge was to identify the significant factors determining the perturbation effects of the Earth 
during the three-burn sequence and select the smallest possible parameter set to adequately rep- 
resent these factors. Three vectors are sufficient to specify the orientation of the most important 
components of the problem: 

hLLO - The angular momentum vector of the initial LLO, together with the assumption of a cir- 
cular orbit at 100 km altitude, defines the geometry of the initial orbit. 

Re - The vector from the Moon to the Earth at the midpoint epoch between TEI-1 and TEI- 3 
(£m)- This epoch is roughly the center of the period of the most pronounced perturbation 
accelerations (near apoapsis). As discussed above, the instantaneous direction of this vector 
determines the direction of the perturbation accelerations. 

Voo - The V-infinity vector of the post-TEI-3 hyperbolic departure trajectory. 

Note that each of these vectors occurs at a different epoch (Iillo at t\. Re at £m » and at £3), 
so they all must be referenced to a common non-rotating frame. 

Fully specifying all of these vectors would require nine parameters, and fully defining their ori- 
entation would require six. However, specifying their orientation relative to one another is possible 
using only three angular parameters. See Figure 7 for an illustration of these angles. Note that these 
three angles begin our parameter list: 

1. a - The relative declination of the vector (i.e., the angle between the initial LLO plane 
and the vector) with positive values when falls on the same side of the plane as 
hLLO- 
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Figure 7 Three-Burn TEI Parameterization Geometry and the Moon-Earth Frozen 
Rotating Frame. The angles a , f3, and 0 define the relative orientation of the three 
vectors: h LLO , Re, and V^. 


2. f3 - The angle between two planes: (1) the plane defined by Re and the vector, and 
(2) the plane defined by Iillo and the vector. The angle is measured in the right hand 
positive sense about the vector from the first plane to the second. 

3. 6 - The angle between the vector and Re. 

The remaining parameters are: 

4. A trEi ~ Time of flight from TEI-1 to TEI-3. This determines the apoapsis distance of the 
trajectory, which will occur in the general vicinity of the midpoint time £m- The perturbation 
effects due to the Earth are expected to be most pronounced near apoapsis, and the apoapsis 
distance will determine (along with the orientation of the orbit and the Earth-Moon distance) 
the magnitude of the perturbation accelerations encountered. 

5. Re - Distance from the Moon to the Earth at this midpoint time. This determines the mag- 
nitude of the perturbation accelerations and, to a lesser degree, the shape of the acceleration 
field. 

6. Voo ~ Magnitude of the vector from the Moon departure hyperbolic trajectory. This is a 
significant factor in determining the orientation of the maneuvers and intermediate orbits, and 
thus partially determines how the perturbations will affect the trajectory. 

7. S - the declination of the vector with respect to the Earth-Moon orbit plane. As the 
Moon-Earth vector approaches, passes, and then moves away from the midpoint epoch, it 
will rotate in the Earth-Moon orbit plane. So this parameter captures information on how the 
Earth perturbation accelerations would vary over the period from TEI-1 to TEI-3. 

Note that parameters 5 and 6 recapture the magnitude of two of the three defining vectors as 
scalars. So in effect, Re and the vector are completely defined in this model, if only in a 
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relative sense. Only the magnitude of Iillo is unnecessary, since this orbit is uniformly specified 
to be a 100 km altitude circular orbit. This particular form of the parameter set has the advantage 
that it retains the parameters (V^o, A trEi, and a) of the simplified version of the TEI Av database 
with no Earth perturbation effects. 

Generation of the TEI Three-Burn Database with Earth Perturbations 


Table 5. 7D TEI Database Parameters 


Parameter 

Range 

Ko 

From 800 to 1450 in 108.33 m/s steps 

Re 

From 356445 to 406706 in 12565.25 km steps 

AtTEI 

From 1 to 2 in 0.5 day steps 

5 

0° 

0 

From 35° to 105° in 14° steps 

a 

From -90° to 90° in 5° steps and ±1°, ±2°, ±88°, ±89° 

P 

From 0° to 360° in 10° steps 


The parameters for the 7D database are given in Table 5. First, epochs were selected which 
represent the different distances (maximum, minimum, and three intermediate distances) from the 
Moon to the Earth. Values for 0 and Voo were selected to correspond to expected ranges of these 
values for the missions being considered here. The parameters a and /3 have the full range of 
possible values. Three flight times, and one declination case was run. For this study, only 6 of the 
parameters are used, since only one 5 case was generated in the database. The maximum value of 5 
is not large (roughly within ±1 1° for the Moon to Earth transit times considered here), and TEI Av 
sensitivity to this parameter was found to be on the order of only 3 m/s. Future database refinement 
may include additional points for this parameter. 

To generate the database, each permutation of the database parameter set must be used to con- 
struct and solve a three-burn optimization problem. This requires the generation of the LLO angular 
momentum vector (Iillo) and Vqo vector for each case. The Vqo vectors generated would not, in 
general, result in trajectories which would return to the Earth. Instead, these vectors are intended 
to bracket all of the real-world cases for which we might need to reference the data. The real- 
world Voo vectors, which would produce trajectories returning to the full range of desired entry 
conditions, should be encompassed by the values of the cases run. 

Since each of the three parameter vectors is defined at a different epoch, they must all be repre- 
sented in the same non-rotating frame. Although technically any non-rotating reference frame could 
be used for the construction of these vectors, it is convenient to use the Moon-Earth Frozen Rotat- 
ing (MEFR) frame (frozen indicating that it is a “snapshot” in time of a rotating frame, to make a 
non-rotating one). Using this frame makes construction of some of the vectors much simpler. The 
MEFR frame is a Moon-centered non-rotating frame defined at the midpoint time of the three-burn 
sequence (/m = (U + £ 3 ) /2), as shown in Figure 7. The x-axis points in the Earth direction, the z- 
axis is along the angular momentum direction of the Earth-Moon system, and the y-axis completes 
the right-handed set. In this frame, the Earth’s position with respect to the Moon is not fixed. Note 
that, in general, tM 7 ^ £2 (i.e., the second maneuver does not occur exactly at the midpoint time). 
Based on the definition of the MEFR frame, the Moon-Earth vector is known for a given epoch to 
be: 

Re = Re^m) [ 1 0 0 ] T (1) 
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where Re(^m) is the magnitude of the Earth-Moon distance at the midpoint time, obtained from 
the ephemeris (in this case, DE 421). The vector right ascension a ^ is computed from 0 and 5 
(using Napier’s rules for right spherical triangles) by: 

aoo = arccos(cos(#)/ cos (5)) (2) 

where right ascension and declination are defined with respect to the MEFR frame. Then, the Voo 


vector (in Cartesian coordinates) is computed via: 

Voo = Foo [ cos(aoo) cos(<5) sin(a 00 ) cos(<5) sin(5) ] T ( 3 ) 

With this information, the following unit vectors can be generated: 

u 3 = Voo/Voo (4) 

U 2 = (Voo xR e )/||Voo xR e || (5) 

Ui = U 2 x U 3 (6) 


These will serve as a framework for construction of the LLO angular momentum vector, which can 
be done as follows: 


hLLO = cos(a) cos(/?)Ui + cos (a) sin(/3)U2 + sin(af)U 3 (7) 

The angular momentum vector can then be converted into the LLO orbit inclination, %llo » and 
longitude of the ascending node, Qllo- 

This completes the inputs to the TEI optimization problem for generation of the TEI Earth per- 
turbation database. This problem is identical to the no-perturbation TEI problem, except that the 
Earth’s gravity is included in the force model (as shown in Figure 1(c)). All the permutations of the 
input parameters resulted in 1,048,950 cases. Each case took an average of 7 seconds to optimize in 
Copernicus. The entire database required about 85 days of CPU time to complete. 

The database was generated by running Copernicus in a series of 7 nested loops (one for each 
of the parameters). The order of the loops was: A tpEE Re* 0, P* ol. In each loop, the 

previous solution is used as the initial guess for the next solution. The a loop is divided into two 
sub-loops (from -90° to 0°, and from 90° to 10°). This aids in convergence, since stepping from high 
to low relative declination magnitude was found to provide better initial guesses than going from 
low to high. Furthermore, the existence of multiple solution families in the three-burn problem in 
the presence of Earth perturbations must be accounted for. In the region of high values of a (i.e., 
high relative declination), there can exist two locally optimal solutions, either of which could be the 
global minimum. Generally, one has a plane change at TEI-2 that is less than 90°, and the other has 
a greater than 90° plane change. Both of these cases are solved for, and the one with the lowest TEI 
Av is selected for the database. 

Each point in the 7D database is then post processed to compute the Av difference (due to the 
Earth gravity perturbations) from the corresponding point in the lower-fidelity 3D database. This 
difference will be called the “Earth perturbation” magnitude (A vep), and is computed by: 

Avep = Av' tei — Avtei ( 8 ) 

Where Av' TEI is the result from the 7D database, and Avtei is the result from the 3D database (for 
the same value of V^, AtpEi , and a). With this, the 7D database becomes an Earth perturbation 


1308 


database, returning the offset from the 3D database result. The reason that this strategy is employed 
(rather than using the 7D database directly to compute the A vtei) is that the 3D database was 
generated at a finer resolution, uses spline interpolation, and better captures the curvatures and 
inflections of the data. The 7D database uses linear interpolation and has a coarser resolution. 
Combining the data in this way provides high resolution on the first-order trends, while the second- 
order Earth perturbation effects are adequately represented by the coarser data. 

Figures 8-11 show TEI Av and Earth perturbation (EP) results from the 7D database. Figures 8 
and 10 show the minimum, mean, and maximum values as functions of each parameter. Figure 8 
shows the sensitivity to TEI time of flight. As expected, while the maximum total TEI Av decreases 
with longer TEI flight times, the maximum EP increases. For the two day flight time, the maximum 
EP is about ±45 m/s. Figure 10 shows the EP variation for the other parameters (for the 48 hr TEI 
transfer time case). As expected, these indicate that the maximum EP also decreases with increasing 
Moon-Earth distance. Interestingly, the a plot shows that for the worst cases ( a near ±90°) the EP 
is always negative, meaning that the 3D database always overestimated the TEI Av for these cases. 
For a near 0°, the EP is always positive, meaning that the 3D database always underestimated the 
TEI Av for these cases. Figure 9 shows representative contours for a specific V^, Re , and 6. Figure 
1 1 shows the variation of the shape of the Earth perturbation contours as 0 varies. 

Querying the Databases 

A schematic of how the databases are used is shown in Figure 12. Multidimensional cubic B- 
spline interpolation 23 is used to query all the databases except the 7D TEI EP database (which uses 
linear interpolation). First, the Voo vector is looked up based on the TEI-3 epoch, the time-of-flight 
(or time-of-flight limits) from TEI-3 to El, and the desired entry conditions (either the El azimuth 
or the entry target polynomials). Once the Voo vector is known, it is used (along with the epoch, 
initial orbit, and TEI flight time) to compute the parameter inputs to the TEI databases. The angles 
a , f3, and 6 can be computed by*: 

a = atan2 [(h L LO • Voo), ||h L LO x Voo||] (9) 

(3 = atan2 [(Re x Voo) x (Iillo x Voo) • Voc/Voc, (Re x Voo) • (Iillo x Voo)] (10) 
0 = atan2[||R E x VooH^Re- V oo)] (11) 

The 3D and 7D TEI databases are used together to compute the TEI Av. The 3D database produces a 
solution A vtei (without Earth perturbations), and the 7D database produces the Earth perturbation 
value Avep- The two are used to compute the three-burn A vtei by: 


A vtei = A vtei ± Avep (12) 

For the expedited return missions, the El azimuth is a free parameter. The databases are queried 
by a global search routine to find the entry azimuth corresponding to the lowest TEI Av. The 
fixed-azimuth Voo database is used in this case. In general, there can be two local minima for 
this problem, which have two different departure Voo vectors at the Moon. These two solutions 
correlate with the two TEI sequence solution families. 


*The atan2 functions use the Fortran convention of a sine, cosine (or Y value, X value) order of inputs. 
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Figure 8 Plot showing how the Earth Perturbation changes for different values of 
A tTEi* As expected, the maximum Earth perturbations occur for longer flight times. 
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(b) Earth Perturbation A vep 


Figure 9 TEI At;(a) and EP Surfaces showing the effect of TEI flight time across 
the parameters a and (3 . All data shown is for the case where = 1125 m/s, Re = 
381576 km, and 0 = 77° 
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Figure 10 Figure shows the maximum, mean, and minimum trends for total TEI 
Av and Earth Perturbation A vep for the specified database parameter. These values 
were calculated across all other parameters. Results shown are for the 48 hr time of 
flight case. 
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Figure 11 Plots showing how the Earth Perturbation changes for different values of 
9. All plots shown are for the case where = 1125 m/s, Re = 368576 km, and A trEi 
= 48 hrs. 
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Figure 12 TEI Database Methodology. This diagram shows how the 6 inputs (blue 
boxes) are used to compute inputs to the 3 databases (yellow boxes) to produce the TEI 
Av result. The computation requires several intermediate variables to be computed. 


RESULTS AND CONCLUSIONS 

The database querying process is very fast (about 5000 queries per second on a desktop PC), 
allowing a very large number of mission cases to be quickly assessed. A scan was conducted for 
each mission case listed in Table 1 of all retrograde departure orbits at 5° resolution in inclination, 
10° resolution in longitude of the ascending node, and 6 hr resolution in time for the 2018-2036 
lunar cycle. Figure 13 shows a snapshot of this data at a single epoch. Figure 14 shows the three- 
burn TEI Av corresponding to the best and worst case departure orbits for each of the mission cases 
for the entire lunar cycle. Table 6 shows a summary of the worst-case epochs over this period for 
each mission case. A comparison is given with the database results and fully-optimized results from 
Copernicus. For each of these cases, the database is within 0.2% of the optimized result. 

The TEI performance database is a powerful tool for mission trades and vehicle performance 
studies and has been used extensively in the Orion project. It allows the determination of the maxi- 
mum TEI Av required to complete a given mission type over the entire lunar cycle and provides a 
statistical view of mission coverage for lower vehicle capability values. It enables a highly-detailed 
overview of the three-burn TEI performance which would not otherwise be available and is being 
used for both vehicle and mission design. 


Table 6. Worst Cases: Comparison of Database with Optimized Results from Copernicus 



Fast 

Expedited A 

Fast 

Expedited B 

Relaxed 

Expedited 

Simplified 

Nominal 

Entry Target 
Nominal 

LLO INC (deg) 

95 

95 

95 

100 

105 

LLO LAN (deg) 

0 

160 

180 

190 

200 

TEI-1 Epoch (TDB) 

11/7/2033 6:00 

1/25/2028 12:00 

3/12/2029 6:00 

4/5/2024 0:00 

4/25/2025 6:00 

Database TEI Av (m/s) 

1264.96 

1233.01 

1156.39 

1239.84 

1236.83 

Copernicus TEI Av (m/s) 

1263.66 

1231.67 

1155.38 

1238.83 

1238.99 

Error (m/s) 

-1.30 

-1.35 

-1.02 

-1.00 

2.16 

% Error 

0.10 

0.11 

0.09 

0.08 

0.17 
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Figure 13 TEI Av contour plots showing retrograde lunar orbit departures for each 
mission case, where TEI-1 occurs on January 1, 2018. 
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Figure 14 Performance Results Over the Lunar Cycle. This plot shows the best 
(lowest A vtei) and worst (highest A vtei) departure orbits at each epoch for each 
mission case. 
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